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Abstract: We present calculations showing that upcoming Cosmic Microwave Back- 
ground (CMB) experiments will have the power to improve on current constraints on 
neutrino masses and provide new limits on neutrino degeneracy parameters. The latter 
could surpass those derived from Big Bang Nucleosynthesis (BBN) and the observationally- 
inferred primordial helium abundance. These conclusions derive from our Monte Carlo 
Markov Chain (MCMC) simulations which incorporate a full BBN nuclear reaction net- 
work. This provides a self-consistent treatment of the helium abundance, the baryon 
number, the three individual neutrino degeneracy parameters and other cosmological pa- 
rameters. Our analysis focuses on the effects of gravitational lensing on CMB constraints 
on neutrino rest mass and degeneracy parameter. We find for the PLANCK experiment 
that total (summed) neutrino mass M u > 0.29 eV could be ruled out at 2a or better. 
Likewise neutrino degeneracy parameters £ Ue > 0.11 and |£„ , | > 0.49 could be detected 
or ruled out at 2a confidence, or better. For POLARBEAR we find that the corresponding 
detectable values are M u > 0.75 eV, £ Ue > 0.62, and |£y / T | > 1.1, while for EPIC we obtain 
M v > 0.20 eV, £„ e > 0.045, and \^ /t \ > 0.29. Our forcast for EPIC demonstrates that 
CMB observations have the potential to set constraints on neutrino degeneracy parameters 
which are better than BBN-derived limits and an order of magnitude better than current 
WMAP-derived limits. 
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1. Introduction 

The CMB is a sensitive probe of basic cosmological parameters such as the spatial curvature 
of the universe and the energy density in baryons, dark matter, and dark energy. Funda- 
mental neutrino properties, such as their masses and the effective number of relativistic 
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degrees of freedom, are already constrained by the CMB. Forecasts for future CMB exper- 
iments, e.g. [1], indicate that neutrino properties will be constrained with unprecedented 
accuracy. These constraints, together with results from next-generation terrestrial exper- 
iments, may enable otherwise unobtainable insights into fundamental neutrino physics. 
These results will be complementary to future terrestrial experiments as well as other cos- 
mological probes (e.g., galaxy surveys [2-5], Lya systems [6,7], joint CMB and galaxy 
surveys [8-11], weak lensing [12-15] and joint CMB and weak lensing [16,17]). 

Most CMB features are imprinted at the epoch of recombination. However, post- 
recombination effects that introduce secondary temperature anisotropy (e.g., lensing of 
the CMB by large scale structure (LSS) and the late integrated Sachs- Wolfe (ISW) effect) 
and polarization (CMB lensing by LSS) can be used to set tighter constraints on certain 
cosmological parameters. Although neutrinos only weakly interact, they have been present 
for the entire history of the universe and can leave their imprint on both the CMB and 
LSS. This allows high-sensitivity and high-resolution CMB experiments to probe neutrino 
properties through the effect of the neutrinos on LSS. 

The impact of neutrinos on the CMB strongly depends on their rest masses. Solar and 
atmospheric neutrino oscillation experiments have shown that at least two neutrino states 
are massive [18]. These neutrino experiments are sensitive to the differences in the squares 
of the neutrino masses, but not to their absolute mass scale (solar: Sm^ = 7.651q' 2 o x 
1CT 5 eV 2 ; atmospheric: [Sm^l = 2.401q;^ x 10~ 3 eV 2 ). In addition, these neutrino 
experiments have not yet determined the sign of Sm^. If it is positive, then the neutrino 
mass states are in the normal hierarchy, with two lighter mass states and one heavier 
mass state; otherwise, if it is negative, then the neutrino mass states are in the inverted 
hierarchy with two heavier mass states and one lighter mass state. To pin down the three 
neutrino masses, a third independent measurement is required, for example, a measurement 
of the total, summed neutrino mass, M v = J2i=i 2 3 Laboratory measurements of the 
neutrino mass-squared differences, imply that at least one neutrino mass must exceed 
0.049 eV. Thus, resolving the neutrino mass hierarchy (i.e., mass ordering of the solar and 
atmospheric mass-squared doublets) may require sensitivity to the total neutrino mass of 
M v < 0.1 eV. For M v > 0.1 eV, the two mass hierarchies are indistinguishable, but if 
the total neutrino mass could be constrained below this level, the inverted mass hierarchy 
would be ruled out. 

Another fundamental issue is how well cosmological probes can constrain the neutral 
lepton number. The lepton number residing in thermal neutrino seas can be characterized 
by neutrino degeneracy parameters, £j = Hi/kBT u (where fii is the neutrino chemical 
potential of the ith species (v e , or v T ), ks is the Boltzmann constant and T u is the 
neutrino temperature) where neutrinos have a Fermi-Dirac distribution. Current CMB 
data does not require the inclusion of neutrino chemical potentials in the cosmological 
model. In standard cosmology, the neutrino degeneracy parameters are assumed to be zero. 
However, there are a number of non-standard mechanisms that could lead to large neutral 
lepton asymmetries [19-22]. Although calculations suggest that these asymmetries may 
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equilibrate in the early universe [23-26], it is interesting to treat the lepton asymmetries 
in the three neutrino flavors independently. 

Extracting neutrino masses from LSS tracers should account for the possibility that 
their chemical potentials do not vanish. A detection of nonvanishing neutral-lepton- 
asymmetry may have far-reaching implications. The current best upper limits on neutrino 
degeneracy parameters, which are invariant under cosmological expansion, are provided by 
a comparison of big bang nucleosynthesis (BBN) calculations with the observed abundance 
of light elements, especially 4 He [27]. Current upper limits on £ from analysis of the CMB 
are of order unity, while upper limits from BBN are on the order £ ~ 0.1. In this work we 
explore how these limits may be tightened by using future cosmological data. 

This paper discusses the constraints on neutrino masses and degeneracy parameters 
that can be obtained from CMB data alone. In particular, we study the experimental 
capacity of PLANCK \ POLARBEAR 2 and EPIC [28] to constrain these parameters. 
We constructed a joint BBN+CMB pipeline which self-consistently solves for the helium 
fraction, Y p , given the other cosmological parameters and allows all three neutrino chemical 
potentials to vary independently of each other. The helium fraction is not an independent 
parameter in our analysis (a similar approach was adopted in [29-31]). Rather, we employ 
a BBN code [32-34] to self-consistently obtain Y p from a given set of other cosmological 
parameters, such as Hq and £ Ue , £,, M and £j, T Here the three neutrino degeneracy pa- 
rameters are treated as phenomenological time-independent parameters, although models 
of time-dependent neutrino chemical-potentials have also been considered in the literature, 
e.g. [22]. Note, however, that neutrino oscillations at the solar mass-squared splitting scale 
can "even-up" the lepton numbers for the different neutrino flavors - a process suggested 
by [25] and shown to work more or less efficiently (to within a factor of ten) by [26], [23] 
and [24] . The ultimate effect of this oscillation-driven process would be to keep all the lep- 
ton numbers all the same and subsequently fixed with time. Y p is an important ingredient 
in the physics of recombination since it determines the Silk damping scale for a fixed baryon 
number. Earlier works discussing the implications of precise CMB observations on helium 
abundance inference are, e.g. [35-37] and more recently [38]. Our analysis also benefits 
from CMB lensing extraction achieved by employing the standard quadratic estimators 
of the lensing potential [39]. This is important in exploring neutrino physics since it has 
been demonstrated that most of the information on neutrino parameters is encapsulated 
in CMB lensing [1,40]. 

This work adds to previous efforts [29-31]) which have attempted to constrain the neu- 
trino degeneracy parameters from CMB or CMB+BBN by including gravitational lensing 
extraction of the CMB allowing the various degeneracy parameters to vary independently. 
Recently, a similar analysis for WMAP5 was carried out which allowed £^ e ^ £„ = £„ T [30]. 
The PLANCK, POLARBEAR, and EPIC experiments have even higher sensitivity and res- 
olution than WMAP. This can facilitate lensing extraction, allowing them to better probe 

x http:/ /www. rssd.esa.int /index. php?project=planck 
2 http: / /bolo. berkeley.edu / polarbear / 
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neutrino parameters. 



The paper is organized as follows. In section 2, we discuss the effects of neutrinos on 
BBN and the growth of structure. Section 3 describes our Monte Carlo Markov Chain 
(MCMC) simulation and the modifications we introduced in CAMB. The degeneracies of 
neutrino mass and helium abundance with other parameters are especially relevant for 
parameter estimations from CMB observations and are therefore extensively discussed in 
section 4. We describe our results in section 5 and conclude in section 6. 



2. Neutrinos and Neutral-Lepton Degeneracy 
2.1 Definitions and Basic Quantities 

Over the history of the universe considered in this work the distribution functions of neu- 
trinos (y) and anti-neutrinos (y) with physical momentum p are 



V t 



U(p;T v ,0 « -r- , (2.1) 

e^ +€ + 1 

where £ = njT v is the degeneracy parameter and T v is the time-dependent neutrino tem- 
perature [41]. From here on, we will use natural units where h = c = &b = 1. The 
degeneracy parameter is a comoving invariant. We assume that at some point in the early 
universe neutrinos and anti-neutrinos were in thermal and chemical equilibrium with the 
photon-baryon plasma and therefore £j, + £p = 0. The cosmic neutrino background (Cz/B) 
temperature is inversely proportional to the cosmological scale factor, a, and is related to 
the (post-recombination) CMB blackbody temperature by T u = (4/11) 1//3 Tcmb- 

It is convenient to write the neutrino energy density and pressure in terms of the 
comoving momentum, q = pa [42] : 

a -4 roc 

Pv + P» = ^2j I 2 <^vV + (aMf [f„(q/a; T„,Z) + f p (q/a; T v ,$] 

4 f oo 2 

P» + P» = ^ Q 2 dq \ =\fv(q/a;T v ,S) + U(q/a;T v ,Q], (2.2) 
ovr^ J y/q? + ( a M) 2 

where M = m u /T U Q, with the Cz/B temperature at the current epoch, T u q 1.95 K. 

The degeneracy parameter is related to the neutral lepton number, 

^^-Ssk^^' <2 ' 3) 

where n u , rip and n 7 are the number densities of neutrinos, anti-neutrinos and photons 
respectively, and £(3) ~ 1.202 is the Riemann zeta function with argument 3. 
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The effective number of relativistic species, N e g, is a ratio between the energy density 
in a given relativistic species and the energy density of the same relativistic species with a 
thermal distribution and zero chemical potential. If a neutrino and anti-neutrino of a given 
flavor both have Fermi-Dirac spectra and have equal and opposite degeneracy parameters 
[e.g., Eq. (|2.1|)], then the effective number of relativistic species contributed by this neutrino 
flavor is 

+"(§)"+"(£)'■ 

The change in the total effective number of relativistic species is often used to describe 
the effects of non-standard neutrino distributions. When different degeneracy parame- 
ters are introduced for each neutrino flavor, the change in the overall effective number of 
relativistic species is 



AiV cff = 



30 /&y 15 

7 \tt 7 Ivr 



(2.5) 



AiV c ff is a useful parameter when the detectable effects of neutrinos on the CMB depend on 
the contribution of these particles to the energy density in radiation. However, upcoming 
CMB experiments will have the sensitivity to probe effects that are dependent on the 
distribution of neutrino energies. In addition, BBN abundances are sensitive to neutrino 
energy distributions. Hence, will be more useful than AA^ in the analysis of upcoming 
CMB experiments. 



2.2 Neutrino Effects on Cosmology 



Neutrinos have a wide range of effects on the evolution of the universe. In the early universe 
they participate in the reactions that determine the neutron-to-proton ratio which, in turn, 
affects the abundances of the light elements produced during BBN. Later, at a redshift of 
z ~ 3200, the energy density in the Ci/B helps determine the epoch of matter-radiation 
equality. At recombination, z ~ 1100, the universe was not purely matter dominated, 
implying that gravitational potential wells had decayed, slightly. This leads to the early 
ISW effect, which boosts the CMB temperature anisotropy angular power spectrum on 
multipole scales associated with the horizon scale, I < 200. Neutrinos play a role in this 
process because the fraction of the total energy density in the form of radiation (which is 
sensitive to neutrino masses and degeneracy parameters, Eq. 2.2) determines the amplitude 
of the ISW effect. This is the only effect of neutrino mass and degeneracy parameter that 
can be probed by WMAP and other moderate angular resolution experiments. Fig. [I] 
shows the calculated CMB power spectrum for various £, along with the data points from 
WMAP5. It is clear that £ > 1 is excluded at la (assuming all other parameters are fixed). 
A global parameter analysis reaches a similar conclusion [30] . 

An aspect highlighted in this paper is that stringent constraints on neutrino mass 
and degeneracy parameters can come from an analysis of CMB lensing. A neutrino that 
is non-relativistic today could have been relativistic at higher redshifts. Non-relativistic 
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neutrinos could be captured into potential wells, while relativistic neutrinos would act as 
hot dark matter (HDM) and would freely stream, resulting in an apparent suppression of 
structure formation during the epochs when the neutrinos are relativistic. Precise mea- 
surements of the LSS power spectrum can be used to place constraints on neutrino masses 
and degeneracy parameters. 

2.2.1 BBN and Light Element Abundances 

BBN occurs at temperatures much higher than the current upper bounds on neutrino 
masses and therefore BBN calculations cannot constrain neutrino masses. However, the 
neutrino degeneracy parameters impact BBN abundance-yields by affecting both the reac- 
tion rates that determine the neutron-to-proton ratio and the expansion-rate of the uni- 
verse, which helps to determine how neutron-to-proton inter-conversion works. The weak 
reactions that set the neutron-to-proton ratio are 

v e + n ^± p + e~ 

v e +p^± n + e + (2.6) 
n ^ p + e" + i> e . 

The rates of these reactions depend on the number density and energy spectrum of v e and 
P e , which in turn depend on the electron neutrino degeneracy parameter, £ Ve [27,43-46]. 
These reaction rates compete with the expansion rate of the universe which is determined 
by the total energy density; the latter also depends on all neutrino degeneracy parameters, 
Eq. (2.2). It is clear, therefore, that BBN distinguishes £ Ve from and £„ r , making for 
a nontrivial interplay between the neutrino degeneracy parameters and the light element 
abundances, particularly 4 He. 

Combined analysis of the CMB (BOOMERANG and DASI experiments), BBN (helium 
and deuterium abundance) and SNIa data yield the following 2a limits [47] 

-0.01 < £„ e < 0.22 

16,^1 < 2.6. (2.7) 

If oscillation between the three neutrino species results in equilibration of the asymmetries 
among the neutrino flavors [23-25], then the more stringent 4 He constraint on £„ e applies 
to all neutrino flavors and BBN considerations suggest [27,45,48] 

16,1 < 0.1. (2.8) 

However, non-standard physics could lead to different degeneracy parameters for the three 
different neutrino flavors. 

BBN determines the abundance of light elements, including the helium fraction, Y p . 
These abundances can be sensitive to the baryon closure fraction, $7^ and the three neutrino 
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degeneracy parameters. In particular, Y p is determined principally by the neutron-to- 
proton ratio at temperatures T ~ 100 keV. This ratio is set by the competition between 
the weak reactions in Eq. (2.6). As a result, Y p depends strongly on the neutrino degeneracy 
parameters. 

Until recently, the helium fraction was usually considered as a free parameter in CMB 
analyses. Recent work [29-31] has attempted to self-consistently include Y p as a non- 
independent parameter in CMB power spectra calculations. It was noted in [31] that 
certain cosmological parameters are significantly biased when Y p is fixed at Y p = 0.24 and 
consistency with BBN is ignored. 

Helium recombination occurs prior to hydrogen recombination. Therefore, for a fixed 
baryon closure fraction, the number density of free electrons at hydrogen recombination 
is a function of the helium abundance. The Silk damping scale is the scale over which 
temperature anisotropy and polarization will be washed out by free-streaming of photons 
between the onset and the end of decoupling. This scale depends on the photon mean free 
path which is inversely proportional to the number density of free electrons. Increasing Y p 
reduces the number density of free electrons at hydrogen recombination, which increases 
the mean free path of the CMB photons. The result would be a suppression of correlations 
on larger angular scales, which would shift Silk damping to lower multipole numbers. 

2.2.2 The Growth of Large Scale Structure 

While weak constraints on neutrino masses can be extracted from the primary CMB power 
spectra, adding probes of structure formation has the potential to significantly tighten 
these bounds. Using CMB lensing rather than resorting to other cosmological probes of 
structure formation is nearly systematic-free, providing high fidelity constraints. 

CMB lensing is a sensitive probe of any cosmological parameter that impacts the 
growth rate of gravitational potential wells. Current CMB data, combined with obser- 
vational data from Type la Supernovae (SNIa) and baryon acoustic oscillations (BAO), 
constrains the total neutrino mass to the sub-eV level [49]. Since the lensed CMB is the 
result of the integrated effect of the lensing of the primary CMB by structure formation, 
and the relevant redshift range for structure formation may overlap with the epoch where 
neutrinos transition from being relativistic to non-relativistic, the CMB can be a powerful 
tracer of neutrino masses and degeneracy parameters. Additional leverage on neutrino free 
streaming comes from e.g., galaxy correlations, Lya forest power spectra [3,50] and weak 
galaxy lensing [16]. 

Tracers of the matter power spectrum, such as CMB lensing, are sensitive to the epoch 
when neutrino momenta were redshifted to a point where they are non-relativistic. This 
is because non-relativistic neutrinos behave as a cold dark matter (CDM) and contribute 
to the growth of structure, while relativistic neutrinos behave as HDM and suppress struc- 
ture on scales below their free streaming scale. Thus, the epoch when neutrinos become 
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non-relativistic is important in discerning the effect of neutrinos on large scale structure. 
Both the neutrino mass and degeneracy parameter determine when neutrinos become non- 
relativistic. 

Probes of the growth of structure in the universe indicate that CDM, rather than 
HDM, is the dominant component of matter. Neutrino masses of 0.2 — 0.3 eV (consistent 
with the upper limits in the current neutrino mass constraints) are mildly relativistic at 
recombination which would result in the slight decay of gravitational potential wells at 
last scattering, leading to a primary ISW effect. For a spatially flat universe, = 0, 
and a fixed dark energy density fraction, Oa, changing the neutrino masses will change 
the amount of HDM at a given redshift at the expense of CDM. This will cause a relative 
suppression of structure formation at high redshifts. In turn, this will be reflected in the 
level of CMB lensing by LSS. Several forecasts for PLANCK, CMBPOL 3 and other CMB 
experiments suggest that constraints on neutrino masses can be improved by a factor of 
three to four [1], provided the experiments have sufficiently high sensitivity and angular 
resolution to allow lensing extraction. As already mentioned, the CMB, galaxy redshift 
surveys, cluster abundances, Lya and other sensitive probes of the growth of structure 
on scales of a few tens of Mpc can be employed to set sub-eV constraints on the total 
neutrino mass. An intriguing question is whether these probes could distinguish between 
the normal and inverted neutrino mass hierarchies. Constraining the total neutrino mass 
below the required ~ 0.1 eV scale is a challenging task in the presence of astrophysical 
foregrounds and other systematics. It was recently shown, e.g. [3,15], that by combining 
several cosmological probes, this (or similar) limit can be achieved. However, it is important 
to be mindful of the assumptions that are made in achieving these limits and to what extent 
the systematics can be controlled. 

The free streaming scale can be estimated as the proper distance traveled by a neutrino 
over the age of the universe. This gives an estimate of the neutrino free streaming scale, 
Afs- This scale is 



where v(t) is the neutrino velocity, which decreases as the universe expands, a(t) is the scale 
factor, (• • •) denotes an average with respect to the neutrino energy distribution function, 



suppressed by neutrino free streaming. This suppression factor will be proportional to Q v , 
the neutrino energy density in closure density units. 

The average free streaming scale for a neutrino species with mass m v and degeneracy 
parameter £ is 






Matter overdensity on scales smaller than Afs will be 



Afs(£,M) 




(2.10) 



F 2 (0 J g=0 Jt =0 eiS + 1 yV + [ a (t)M] 2 a(t) 



3 http:/ /cmbpol. uchicago.edu/ 
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where M = m u /T V Q and is the Fermi integral of order two, 

The free streaming scale dependence on neutrino masses and degeneracy parameter is illus- 
trated in Figure ||. Note that the free streaming scale increases with decreasing mass and 
increasing degeneracy parameter. Care should be taken when simultaneously discussing 
neutrino degeneracy parameters (which are related to flavor states) and neutrino masses 
(which are related to mass states) [41]. 

Figure || illustrates how neutrino mass and degeneracy parameters affect the transfer 
function. The transfer function represents the effect of all physical processes that cause 
the primordial power spectrum to evolve into the matter power spectrum at latter epochs. 
The relation between the power spectrum P m {k;z) and the transfer function T(k;z) can 
be written as 

P m (k;z) = A s k n °T 2 (k;z), (2.12) 

where n s and A s are the tilt and normalization of the primordial power spectrum. The 
suppression of the matter power spectrum on scales smaller than the neutrino free streaming 
scale is related to observable quantities. This could be obtained from galaxy surveys 
or inferred from the CMB angular power spectrum by deconvolving the lensing power 
spectrum. The change in the matter power spectrum resulting from neutrino free streaming 
is [51] 

^Pm(k) ^ (2 13) 

Here, P m (k) = P m {k\z = 0). The effect of non- vanishing neutrino mass is shown on the 
left side of Figure ||. As the neutrino mass increases, the free streaming scale decreases, 
leading to suppression of the transfer function at larger wavenumbers (note that all curves 
are normalized to the case of m v = at low k). The suppression of the transfer function at 
these large wavenumbers is more pronounced for larger neutrino masses because in this case 
neutrinos constitute a larger fraction of the dark matter. The effect of non-zero neutrino 
degeneracy parameter is shown on the right side of Figure [3|. As the degeneracy parameter 
increases, the free streaming scale increases, leading to suppression of the transfer function 
at progressively smaller wavenumbers. 

Figure |] illustrates the effect of non-zero degeneracy parameters on the CMB temper- 
ature, polarization, and deflection angle power spectra (the latter is essentially a measure 
of the rms lensing deflection angle of the LSS, relevant to CMB lensing). The most signifi- 
cant differences are at large multipoles, to which current CMB experiments are blind, but 
PLANCK, POLARBEAR and EPIC will be sensitive. Although degeneracy parameters 
= 3 are already ruled out by BBN and CMB data, we show these cases for illustrative 
purposes. From the plots of Cf T and Cf E we can see the effect of neutrino degeneracy 
parameters on scales from the acoustic horizon at recombination down to Silk damping 
scales. The power spectrum Cf d for lensing deflection angle, d, is suppressed in the pres- 
ence of nonvanishing at high I. This reflects the relative suppression in the transfer 
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function at large wavenumbers. This effect, leads to suppression of the lensing-induced 
B-mode polarization that results from E-B conversion via CMB lensing by LSS. Therefore, 
the presence of a large £„ can be constrained with high sensitivity CMB experiments which 
are capable of discerning small variations in the weak B-mode polarization. 

3. CMB Code and Monte Carlo Simulation 

As in [29,31,42], we modified the Boltzmann code CAMB [52] by replacing the neutrino 
distribution function (here q=p/T) 

with 

everywhere, including in the expressions for energy density and pressure, as well as in the 
Liouville equation for neutrino density perturbations. We allowed the individual neutrino 
flavors to have three different degeneracy parameters £„ e , and £ Vt . 

Neutrino masses are subject to experimental constraints from terrestrial neutrino ex- 
periments. The neutrino masses and the degeneracy parameters are degenerate with various 
other cosmological parameters, presenting a challenge in attempting to determine these pa- 
rameters using the CMB. For a given value of both m v and are degenerate; for a 
fixed Qv, increasing must be compensated by decreasing neutrino masses. To avoid this 
degeneracy in interpretation of our simulation results we consider mi, m-2, ms, £ Ue , £ u and 
£j, T as our basic parameters (in addition to the standard cosmological parameters). The 
three neutrino masses are constrained by the measured mass squared differences. In this 
work, we used conservative gaussian priors for these differences: Sm^i = 8.0±0.6x 10~ 5 eV 2 
and Sm^i = 2.4 ± 0.6 x 10~ 3 eV 2 , which are consistent with the experimental values. In 
practice, these and projected future laboratory improvements in neutrino mass-squared 
uncertainties, have little effect on our analysis. The dominant uncertainties come from 
CMB data uncertainties. 

CMB data is encapsulated in the angular power spectra down to scales determined 
by the angular resolution of the specific experiment and its instrumental noise level as 
compared to the CMB signal. The instrumental noise Ni in measuring the angular power 
spectrum for multipole /, for the autocorrelation of the temperature and polarization of 
the E- and B-modes are related for bolometric radiometers by 2Nf T = Nf E = N^ B . The 
instrumental noise is uncorrelated between T, E, and B and increases exponentially with 
multipole number, 

JV/J = ,5 ab (# a A a ) 2 exp[Z(Z + l)0 2 /81n2], (3.3) 
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where a and b are either T, E, or B. Here, the noise at the frequency band centered at 
v is a function of the corresponding beamwidth, 8 a , and the noise per pixel in equivalent 
temperature units, A . To obtain the effective noise power contributed by all frequency 
bands in the experiment one adds them as if they were uncorrelated gaussians 

jvt = ew 1 ]" 1 - ( 3 - 4 ) 

V 

We simulated parameter extraction from PLANCK, POLARBEAR and EPIC. The 
sensitivity and resolution we considered for these experiments are given in Table [l]. The 
CMB power spectra C TT , C TE , C EE and C BB , together with the power spectrum of the 
deflection angle C dd , and its cross-correlation with the temperature anisotropy, C Td are 
calculated by CAMB for a given fiducial cosmological model. When the Cf d power spectra 
are calculated with the parameter accuracy_level=l, there are very large oscillations at 
I ~ 200. These go away when the accuracy level is increased. In order to fix this in our 
simulation, the Cf d at I > 200 in the simulated data are replaced with values calculated 
when setting accuracy _level=5 in CAMB. When calculating the likelihood in CAMB, the 
proposed Cf are set to the same value as the "experimental" Cf d . We found that doing 
so does not affect the parameter uncertainties extracted from the MCMC simulation. All 
power spectra are assumed gaussian and are taken to be unlensed following the conclusion 
of [1]. The noise in lensing reconstruction, A 7 , , is a function of the observed power spectra 
(all four lensed C\ with instrumental noise, Eq.(3.3), included) and the unlensed power 
spectra [40] (without lensing, as obtained from CAMB for a fiducial cosmological model). 
In calculating N dd we employ the publically available code [53] which makes use of the 
quadratic estimators [39]. 

4. M v and and Their Degeneracies with Other Parameters 

Neutrino masses and chemical potentials are both degenerate with each other and with 
other cosmological parameters. We now discuss the main degeneracies of neutrino pa- 
rameters with other cosmological parameters. In particular, we discuss the degeneracy of 
neutrino masses in section 4.1 and chemical potential in section 4.2. 

In Figures |5| - [ll], we examine the degeneracies between different cosmological parame- 
ters and parameters that are affected by neutrino physics, in particular M u , £ u and Y p . In 
these figures, we assume the normal mass hierarchy with m\ = 0.01 eV, which corresponds 
to M u = 0.073 eV when the priors on the neutrino mass squared differences mentioned in 
the previous section are imposed, and £ v = when calculating the theoretically expected 
CMB power spectra. All other cosmological parameters are set to the best fit values of 
WMAP [49]. As a result, these figures illustrate the possible constraints that could be 
placed on the neutrino parameters if the actual values of these parameters are too small 
to create an effect on the CMB that could be discernible in the upcoming CMB experi- 
ments. One important point to keep in mind is that the neutrino mass-squared differences 
measured in the laboratory excludes any M v < 0.058 eV. 
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4.1 Neutrino Mass Degeneracy 



Cosmological probes of neutrino masses are sensitive to the kinematics of individual neu- 
trinos. Since the gravitational interaction is flavor-blind, it does not distinguish between 
neutrino species. However, for a fixed total mass it does depend on how this mass is 
distributed between the three species [42] . 

4.1.1 Degeneracy with w 

As mentioned above, the suppression of the matter power spectrum in the presence of 
massless neutrinos on scales much smaller than the neutrino free-streaming scale is 
AP m (k)/P m (k) ~ — 8fJ„/fi m . Thus, increasing Vt u as a result of increasing the neutrino 
mass can be compensated by increasing f2 m . However, to keep the universe spatially flat, 
the closure fraction of the dark energy must be lowered - this is achieved by forcing the 
dark energy equation of state parameter, w, to be more negative. Therefore, increasing m v 
is degenerate with lowering w as illustrated in Figure ||. One way to avoid this degeneracy 
is to ignore cosmological information from scales smaller than the neutrino damping scale 
(which comes at the cost of significantly weakening the power of the CMB as a diagnostic 
tool of neutrino properties). Another possibility to avoid the M v — w degeneracy is to 
employ supplementary measures of distance, e.g., BAO or SNIa [49]. 

4.1.2 Degeneracy with erg 

The fluctuation in the matter density on 8/i -1 Mpc scales is 



where W(kR) is a window function, h = iTo/(100 km/s/Mpc), and R = 8/i _1 Mpc. There- 
fore, as is a function of both A s and n s , the normalization and tilt of the power spectrum, 
respectively. It is also a function of neutrino masses and degeneracy parameters, as well 
as any other cosmological parameters which may affect structure formation and the evo- 
lution of LSS on scales smaller than a few Mpc (Eq. 2.12). Since erg represents the mass 
fluctuation on ~ 10 Mpc scales, and is therefore subject to neutrino free-streaming we can 
expect a slight M v — og degeneracy, at least for CMB experiments which are sensitive to 
angular scales that correspond to neutrino free streaming scales. However, this degeneracy 
is very weak in practice as can be seen from Figure 6; no such degeneracy is expected to 
be observed in the PLANCK data. 

4.1.3 Degeneracy with Hq 

Previous studies have shown an anti-correlation between neutrino mass and the Hubble 
constant. The anti-correlation results from the fact that while all three neutrino mass 
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states are at least mildly-relativistic at recombination (consistent with the WMAP con- 
straints on neutrino mass [49]), at least two of these mass states are non-relativistic today 
(consistent with the 5m? values from neutrino experiments). As a consequence, the neutri- 
nos contribute to Cl m today, but contributed to Q r (the closure fraction in radiation energy 
density) at recombination. Since neutrinos with larger rest masses will constitute a larger 
fraction of the CDM at the current epoch, larger neutrino masses imply a larger O r , relative 
to Q m , on the surface of last scattering which, in turn, gives an enhanced ISW effect. Most 
of the extra power due to this effect is on scales somewhat larger than the horizon (the first 
acoustic peak), effectively extending the first acoustic peak to larger scales. This effect can 
be mimicked by lowering Hq, because a lower Hq implies a larger horizon at decoupling. 

Figure [/] shows the degeneracies between M v and Hq for PLANCK, POLARBEAR, 
and EPIC. In these plots, the M v — Hq degeneracy described above is not evident. This 
results from the fact that these high resolution experiments will constrain neutrino masses 
primarily from lensing information, instead of through the ISW effect at low multipoles. For 
these high resolution experiments, the neutrino free streaming length is the more relevant 
quantity that relates to observables. This is an example of how the higher resolution and 
sensitivities of upcoming CMB experiments open windows to new effects that may lift 
parameter degeneracies. 

4.2 Neutrino Chemical Potential Degeneracy 

Including nonzero neutrino degeneracy parameters in the analysis introduces new param- 
eter degeneracies. In Figure ||, the degeneracy between £ v and M„ is shown. A naive 
interpretation is that for a given neutrino free streaming length, increasing £ w must be 
compensated by increasing the neutrino mass. However, the range of allowed neutrino 
masses and chemical potentials does not allow such a parameter degeneracy in the neu- 
trino free streaming scale (see also Fig. |2|). The mild degeneracy shown here comes from 
the physics at recombination through the decay of potential wells and the ISW effect. 
Nonvanishing neutrino degeneracy parameters increase the energy density in neutrinos for 
fixed neutrino masses. Increasing the neutrino energy density must be compensated by 
increasing the density of CDM in order to keep AP m /P m « — 8Q U /Q m unchanged. This 
degeneracy can be seen in Figure |9[ 

4.3 Helium Fraction Degeneracy 

The helium fraction affects the physics of recombination primarily by changing the Silk 
damping scale. The baryon closure fraction, O5, is obtained to high precision from the 
amplitudes of the acoustic peaks of the CMB. For a given Qf,, more helium implies less 
hydrogen and fewer free electrons on the surface of last scattering. This causes a larger 
photon mean free path, which damps CMB temperature anisotropy on larger angular scales. 
This effect can be mimicked by either reducing the normalization of the primordial power 
spectrum or increasing its tilt. Both the Y p — A s and the Y p — n s planes are shown in 
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Figure 10. There is a significant difference in the Y p axes between the £ v = and the 
£ v 7^ cases; BBN data and the current precision on cosmological parameters tightly 
constrain Y p , but allowing £ v to be nonzero affects the BBN-calculated yield for Y p and 
allows these degeneracies to manifest themselves in the analysis. The dilution of the free 
electron density at the epoch of last scattering also can be compensated by increasing l^/i 2 
as can be seen in Figure O. 



5. Results 



We adopt a 14 parameter cosmological model with priors on neutrino masses taken from 
neutrino oscillation data. Our model is consistent with the concordance cosmological model 
[49]. With this model, a BBN+CMB MCMC analysis demonstrates that PLANCK and 
POLARBEAR will be able to measure a total neutrino mass of 0.29 eV (PLANCK) and 
0.75 eV (POLARBEAR) at the 95% confidence level. In addition, neutrino degeneracy 
parameters can be constrained to be smaller than 0.11 (£ e ) and 0.49 (£„, £ T ) for PLANCK 
and 0.62 (£ e ) and 1.1 £ T ) for POLARBEAR. The former constraint on £ Ve is already 
better than the corresponding BBN one, Eq. (2.7). 

It is interesting to examine the sensitivities of upcoming CMB experiments to the 
neutral lepton asymmetries. If the neutrino asymmetries equilibrated through neutrino os- 
cillations prior to the BBN epoch, then the possible constraints on the neutrino degeneracy 
parameters become very strong, £ Ve = £ u = ^ Vt < 0.06 [54]. 

Our analysis with the extended parameter space yields weaker constraints on neutrino 
masses. For example, in the minimal model (all are set to 0) the PLANCK 2a upper 
limit on M v is 0.27 eV. When the degeneracy parameters are turned on, the corresponding 
upper limit on M u rises to 0.29 eV. The la confidence range on the electron neutrino 
degeneracy parameters are -0.0314 < £„ e < 0.108 for PLANCK and -0.017 < £„ e < 0.33 
for POLARBEAR. The reason for the skewness of the distribution towards positive £ Ue 
values results from the fact that v e determines the reaction rates of the processes described 
in Eq. (2.6) and also the expansion rate. In contradistinction, £ v and £„ T affect only the 
expansion rate. 

Perhaps the ultimate CMB experiment to address B-mode related issues is the mission 
concept, EPIC. We find that EPIC will be able to set an upper limit on the total neutrino 
mass of ~ 0.20 eV at 2a confidence. In addition, we find that EPIC data alone will have 
sufficient sensitivity to constrain the neutrino degeneracy parameters to a level which can 
compete with the current constraints on degeneracy parameters derived from the primordial 
abundance of light elements. We derive the following limits on the degeneracy parameters: 
£i/ e < 0.045 and £„ T < 0.29 at 2a. These are better than current BBN constraints, 
even if equilibration of the degeneracy parameters is assumed. These results are achievable 
without resorting to assumptions about flavor mixing in the early universe. EPIC is capable 
of such an improvement in sensitivity to the degeneracy parameters because of its very high 
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sensitivity and angular resolution which allow for the precise measurement of the B-mode 
polarization required for lensing extraction of the CMB. 

Finally, a cautionary note on the efficacy of our analysis. It is reasonable to ask 
whether future CMB data will indeed warrant considering a 13- or 14-parameter model. 
This question is reasonable even in the context of idealized analysis presented here because 
CMB data is at least limited by cosmic variance and instrumental noise. In the real world it 
will also be limited by astrophysical foregrounds as well as systematics. As in [55], adding 
cosmological parameters is in general expected to improve the fit of data to the theoretical 
model. Defining a generalized x 2 such as \ 2 = — 21n(L) + 2p (where L is a likelihood and p 
is the number of cosmological parameters in a given model) and exploring if it improves is a 
useful test for such models [55]. However, this requires real data, i.e. sky-maps in the CMB 
case. Our analysis employed a mock power spectrum which in principle we could use to 
generate multiple sky realizations, each one yielding a different numerical value of x 2 . One 
can then statistically determine what fraction of these actually improve when we extend 
the model from 11 to 13 or 14 parameters. However, this procedure is time consuming 
and may not be necessary at this point. When the real CMB data considered in this work 
is available it will be straightforward to determine whether or not our generalized model 
gives a better fit to the data. 

6. Conclusion 

Within a decade the CMB has transformed from being a cosmological probe of the basic 
cosmological parameters to a probe of physics beyond the standard model. CMB experi- 
ments have set interesting limits on the energy scale of inflation as well as on exotic physics 
such as topological defects from phase transitions in the early universe and cosmological 
birefringence. Also, important constraints on neutrino masses have already been obtained. 
Although WMAP has constrained neutrino masses to the sub-eV level, it is an exciting 
possibility that precise CMB measurements could place stringent constraints on neutrino 
masses and neutral lepton asymmetries. Although the Cz/B neutrinos cannot be directly 
detected, they can be indirectly detected through their dynamics (through their effect on 
the expansion rate) and kinematics (via the damping of LSS by neutrino free streaming). 
A convincing detection of the CuB would be a monumental discovery in the history of 
cosmology. 

In this paper we explored the effects of neutrino mass and nonzero neutrino degeneracy 
parameters on the CMB. Changing £ u leads to a different value of Y p from BBN, and Y p 
affects the density of free electrons at recombination. Y p was calculated self-consistently in 
our analysis by a BBN code with a given set of cosmological and neutrino parameters. No 
priors on Y p were included or used in the analysis. 

Our analysis is conservative in that it allows the three neutrino degeneracy param- 
eters to be independent parameters in the analysis. Typically, the neutrino degeneracy 
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parameters are assumed to be equal to each other, which would be expected if the neutrino 
asymmetries equilibrated in the early universe. Other works have at least set = £„ T , 
since the physics of the early universe is insensitive to the difference between these param- 
eters. However, neutrino free streaming lengths are sensitive to the absolute value of each 
neutrino degeneracy parameter, so we treated each degeneracy parameter independently. 
While the addition of data from other cosmological probes of distance scales or LSS could 
be included to break parameter degeneracies, we did not include them so that we could 
isolate the probative powers of the CMB alone. 

Upcoming CMB experiments such as PLANCK, POLARBEAR, and perhaps also 
EPIC will have the capability to either detect the neutrino masses and degeneracy param- 
eters or place much more stringent bounds on these parameters as compared to current 
constraints from WMAP. These breakthroughs in the power of the CMB to detect neu- 
trino parameters is the direct result of the improved resolution and sensitivity of upcoming 
CMB experiments, allowing CMB lensing extraction to provide an ultra-sensitive handle 
on neutrino masses and degeneracy parameters. 
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Experiment / s k y 
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3.2 
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13.3 




850 


1.0 
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1047 



Table 1: Sensitivity parameters of the CMB experiments considered in this work: / s k y is the 
observed fraction of the sky, v is the center frequency of the channels in GHz, 9b is the full width 
at half maximum in arc- minutes, At is the temperature sensitivity per pixel in [iK and = As 
is the polarization sensitivity. 
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Figure 1: The calculated CMB temperature anisotropy power spectrum for £„ = (fiducial model), 
0.5, 1.0, and 2.0. The WMAP5 data points are included for reference. 
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Figure 2: Neutrino free streaming scale: The figure on the left is the neutrino free streaming scale 
as a function of neutrino mass with £„ = 0. The figure on the right is a contour plot of constant 
free streaming scale in the m v -£ plane; the contours from right to left correspond to Afs = 0.8, 1.2, 
1.6, and 2.0 Gpc/h. 
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Figure 3: Susceptibility of the transfer function to M v (left) and (right). The values used for 
£ are 0.1 (blue), 0.5 (cyan) and 1.0 (yellow). 





Figure 4: CMB power spectra response to changing £„: Cf T (top- left), Cf E (top-right), Cf d 
(bottom). In all three plots the black curves correspond to the fiducial model (£„ e = = £„ T = 0), 
the red curves correspond to a non-zero £„ e model (£„ e = 3, = £„ T = 0), and the blue curves 
correspond to a non-zero £„ T model (£ W(! = 0, = £„ T = 3). 
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Figure 5: The M v -w degeneracy: Shown are the results from PLANCK £„ = (top left), ^ 
(top right), POLARBEAR ^ (bottom left) and EPIC ^ (bottom right) simulations. In 
this plot and each successive plot, the contours correspond to the 1- and 2-a regions. . The fiducial 
cosmological model is WMAP best-fit data and the neutrino masses mi and 7713 subject to neutrino 
oscillation results with mi assumed O.OleV. 
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Figure 6: The M v -cr% degeneracy: Shown are the results for PLANCK in the case £„ — (left) 
and £„ 7^ (right). The fiducial cosmological model is as described in Fig. 5. 
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Figure 7: The M^-Ho degeneracy: Shown are PLANCK (£ = and £ ^ cases on the top-left 
and top-right, respectively), POLARBEAR (left bottom) and EPIC (right bottom) results. For 
POLARBEAR and EPIC we show the generalized cosmological model with £ v ^ 0. The fiducial 
cosmological model is as described in Fig. 5. 
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Figure 8: The M„-£„ degeneracy for PLANCK: The fiducial cosmological model is as described 
in Fig. 5. 
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Figure 9: The - n c h 2 degeneracy: Shown are results for PLANCK (top left), POLARBEAR 
(top right) and EPIC (bottom). The fiducial cosmological model is as described in Fig. 5. 
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Figure 10: The degeneracy of Y p with the normalization and tilt of the primordial power spectrum. 
The results for PLANCK with = are on the left and with ^ are on the right. The top 
two plots depict the Y p -A s degeneracy and bottom plots show the Y p -n s degeneracy. The fiducial 
cosmological model is as described in Fig. 5. 
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Figure 11: The Y p — f2&/i 2 degeneracy: Shown are results for PLANCK with = (top left), 
PLANCK with ^ (top right), and POLARBEAR with ^ (bottom). The fiducial cosmo- 
logical model is as described in Fig. 5. 
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